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Abstract 



The problem of particle oscillation is considered in a pedagogical and compre- 
hensive way. Examples from K, B and neutrino physics are given. Conceptual 
difficulties of the traditional approach to particle oscillation are discussed. It 
is shown how the probability current density and the wave packet treatments 
of particle oscillations resolve some problems. It is also shown that only full 
field theoretical approach is free from conceptual difficulties. The possibility 
of oscillation of particles produced together with kaons or neutrinos is con- 
sidered in full wave packet quantum mechanics language. Precise definition 
of the oscillation of particles which recoil against mixed states is given. The 
general amplitude which describes the oscillation of two particles in the final 
states is found. Using this EPR-type amplitude the problem of oscillation of 
particles recoiling against kaons or neutrinos is resolved. The relativistic EPR 
correlations on distances of the order of coherence lengths are considered. 
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I. INTRODUCTION 



The subject is not new. The problem is known since 1955 when Gell-Mann and Pais [1] 
predicted the existence of two neutral kaons. Earlier, in 1953 in a scheme for classifying the 
various newly-found particles, Gell-Mann represented the neutral kaon K° and its antipar- 



ticle K as two distinct particles. The decay of both particles into 7r + 7r~ was observed. If so, 
how do we know which particle has originated it: the K° or the K ? The problem has been 
solved by realizing that what we observe is the mixture of two states, K° and K : 
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where e is a small, complex, later measured parameter responsible for CP symmetry breaking 
[2]. In this way first time the interference between states of slightly different masses has 
appeared in quantum mechanics. Inspired by the work of Gell-Mann and Pais, Bruno 
Pontecorvo turned to consider the possibility of quantum mechanical mixing in another 
neutral particle - the neutrino. In 1957 he first suggested that a neutrino may oscillate into 
its antipartner [3] . Oscillation among the different kinds of neutrinos was then proposed by 
Maki, Nakagawa and Sakata in 1962 [4] and later by many others [5]. 

The neutral K° — K boson system is not the only one where the quantum mechanical 
mass mixing can be considered. We can expect to observe the same phenomena in other 
neutral boson systems: D° — D and B° — B . Generally, flavour oscillations of particles 
can occur when states produced and detected in a given experiment, are superpositions 
of two or more eigenstates with different masses. The oscillation of K and B meson has 
been observed experimentally in 1961 [6] and later [7] and has been used to place stringent 
constraints on physics beyond the Standard Model. If neutrinos are massive and oscillate 
it is possible to resolve the well-known solar neutrino problem [8]. There are also first 
experiments in which the neutrino oscillations are observed [9]. 

The flavour oscillation of particles is a very fascinating demonstration of quantum me- 
chanics in the macroscopic world. It has served as a model for many interesting systems and 
problems. Various aspects of quantum mechanics, as for example coherence, decoherence, 
wave packets, measurements, similarity and differences between pure and mixed states, wave 
function collapse, EPR "paradox" are in action. On the other hand particle mixing is the 
place where fundamental symmetries and properties of fundamental interactions are stud- 
ied. Discovering of the CP symmetry violation and the measurement of differences between 
neutral mesons masses are connected with K, B bosons mixing. Neutrino oscillations have 
a chance to be the first place where problem of neutrino masses can be resolved. 

In this review we will concentrate only on the quantum mechanical description of particle 
oscillations. Problems connected with testing of the fundamental interactions will not be 
discussed. 

First of all we should mention that interference between states with different masses is not 
allowed in non-relativistic quantum mechanics. The Galilean invariance forbids a coherent 
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superposition of such states (the so called Bargman superselection rule [10]). Beyond the 
non-relativistic limit such restrictions do not hold (which clearly follows from experiment). 
It means that all of our considerations should be done in relativistic quantum mechanics 
(nevertheless non-relativistic approximations are possible). 

First, we would like to describe briefly the traditional approach to the particle oscilla- 
tion problem. This treatment is simple and elegant but immediately raises a number of 
conceptual questions. We specify more of them (Chapter 2). Next we show the wave packet 
treatment, where some of the problems disappear (Chapter 3). The current density approach 
which is closely connected with the experimental setting, is described in Chapter 4. The 
problem of constructing the probability current density for a particle with undetermined 
mass is also considered there. 

Next, in Chapter 5, we give some remarks on the field theoretical approach to particle 
oscillations. Usually (as in the case of neutrino oscillations) the oscillating particle is not 
directly observed. Only particles accompanying neutrinos, hadrons and charged leptons 
created in the decay are observed. The proper approach should take all these circumstances 
into account. The creation of the neutrino in the source, its propagation to the detector 
and the detection process are treated in the framework of quantum field theory as one large 
Feynman diagram. 

In Chapter 6 we discuss the controversial problem of the oscillation of particles recoiling 
against kaons or neutrinos from the production process. A detailed approach using wave 
packets explains the problem of four-momentum nonconservation raised in the literature. 

In Chapter 7 we discuss the modern example of an Einstein-Podolsky-Rosen correlation 
in K° — K and B° — B systems. The amplitude approach does not entail the somewhat 
mysterious "collapse of the wave function" which is usually invoked to describe the EPR 
effects. 

Finally in Chapter 8 we summarize our main conclusions. 



II. PROBLEMS CONNECTED WITH THE TRADITIONAL APPROACH TO 

THE PARTICLE OSCILLATION. 

The usual description of kaon mixing phenomena can be found in many textbooks [11]. 
Suppose, that we produce K° at t = by the reaction 



n-p -> K° A°. 

From (1) and (2) the K° state at t = is 
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After time t, as |fCs)and \Kl) states are definite mass eigenstates, we have 
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where mug) an d ^l(S) are masses and inverse mean lifetimes respectively of the long (short- 
lived component of K. 

The K° (K° ) fraction of the beam after time t is just 
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(t) = (K u (W )\K» (t) 
1 
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(K° )\K° 
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where Am = tjil — mg. From Eq.(6) we can see that the fraction of K° (K° ) becomes 
smaller (because of decay) and changes with time with frequency uo = • 

Neutrino oscillations are described in a very similar way [12]. Let us assume that at 
t = neutrino with flavour a was born with momentum p perfectly defined (as for example 
neutrino in the pion decay n + — > /i + z^t). At this time the neutrino state is described by 



l^(o)) = E^|a> 



(7) 



where states \a) are energy- momentum eigenstates for neutrinos with mass m a and U aa are 

elements of a flavour-mass mixing matrix. 

Then 



H\a)=E a \a) 



(8) 



where E a = yj p 2 + m 2 with the same momentum p for each neutrino. After time t the state 
will evolve into 



|*a(0)> - !*«(*)> = e~ iHt |tf a (0)> = E^e"* ( \a) . 



(9) 



Then the probability that a neutrino born at t — with flavour a at time t has flavour (3 is 
given by 



p a ^(t) = \(yp(o)\y a (t))t 



E^V 



-iE a t 



a=l 



(10) 



where n is the number of interfering light neutrinos[|. 

Now usually relativistic approximations are made. As for real, light neutrinos p 3> m a , 
we have 

(i) £ a =P+f , 
and 

(ii) a neutrino born in x = 0, at time t will be approximately at position x m t. 



1 To prove Eq.(10) we have to assume that the scalar product of two eigenmomentum states 
(b | a) = 5 a b- This means that we must introduce some normalization volume and momentum and 
energy are quantized. 
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Then from (10) we can find that the probability for our neutrino, born with flavour a, 
to have new flavour (5 after traveling a distance x, is 



Pa^(x) = £ \U Pa \ 2 \U aa \ 2 + 2Y,\u; a U aa Up b U* ab \ cos (2tt-^ - Vm ) , (11) 

a=l a>b ^ Lah ' 



(12) 



a>b 

where L ab , known as oscillation length between u a and v b is defined by 

m a ~ m b 

and 

¥ab-,*p = arg (ui a U aa U pb U* ab ) , (13) 

are phases responsible for CP violation. 

From (11) it follows that the oscillation will disappear (the P a ^p{x) does not depend on 
x) if (i) all neutrino masses are equal m a = m b and/or (ii) only diagonal elements of the 
mixing matrix U ab do not vanish. 

The presented arguments seem to be clear and elegant but they are wrong. Many concep- 
tual questions arise when we look at the presentation shown above. A complete treatment 
of particle oscillation must address the following additional issues. 

(1) A necessary condition for particle oscillation to occur is that particle source and detector 
are localized within the region Ax much smaller than the oscillation length \L ab \. 

\L ab \ > Ax. (14) 



(2) From Eqs. (4) and (7) we see that different mass eigenstates are produced and detected 
coherently. This is possible only if the momentum (p) and energy (E) of the oscillating 
particle are spread in such a way that the error in m 2 measurements given by 



Am 2 



(2E) 2 {AEf + (2p) 2 (Apf 



1/2 



is larger than \m 2 a — m\\ = |Am^ 5 |, 



Am > I Am 



ab I 



If this condition is not satisfied and \m 2 a — m b \ > Am 2 , then also 

> 2pAp. 



(15) 



(16) 



(17) 



But from the uncertainty relation Ax > and Eq.(17) gives AX > n-^r = Which 



Ami 



is in contradiction with Eq.(14). 

From both conditions described above we see that the oscillating particle state cannot 
be described by a plane wave with definite momentum [13] and the wave packet approach 
must be constructed. 
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(3) The energy and momentum conservation in processes in which oscillating particles are 
created (e.g. 7r~p — > AK° or 7r + — > /U + z/^) implies that different mass eigenstate components 
have different energy and momentum [14] . Approaches where all oscillating particles have the 
same momentum and different energies [11,12], or the same energies and different momenta 
[15, 16, 17] are conceptually not correct. 

(4) In the traditional approach to find the oscillation probability we calculate the wave 
functions overlap (compare Eqs. (6) and (10)). This procedure gives the probability which 
depends on time. In the real experiment the distance between the source and the detector is 
known (not the moment in which the measurement is done). To transform P(t) into P (x) 
the classical formula x — vt is invoked. However to find the probability that the beam 
of particles produced at x— will reach a physical detector at a distance 

density j (x,t) should be integrated over the surface of the detector and over the time of 
observation [18] 

,h <t <t 2 ) = dt I dS-~j(x,t). (18) 

' Jti JdA 

In such an approach there is no problem of how to change t into x. 

(5) In case of neutrinos additional conceptual problems arise because neutrinos are not 
"seen" directly. The only things which can be "seen" are hadrons or/and charged leptons in 
points where neutrinos are produced and detected. So in a realistic description, the external 
(initial and final) particles should be described by wave packets, and the mass-eigenstate 
neutrinos should propagate from the production region to a detector [19, 20]. 

All points which have been mentioned above are not only purely academic. We do not 
try to derive in a more precise way something which is known from the beginning. We will 
see that the more precise approach to the particles' oscillation phenomenon gives us new 
predictions and elucidates mysteries in many points. On the other hand we will see quantum 
mechanics in action on macroscopic distances. 



x the current 




III. THE WAVE PACKET TREATMENT OF PARTICLE OSCILLATION. 



The wave packet approach to neutrino oscillation was first proposed by B. Kayser [13] 
and later considered in more detail in [21, 22]. Nowadays the neutral boson oscillation is 
also treated in the same way. We will present the formalism for neutrinos, but everything 
can be repeated also for bosons. 

In Eq.(7) states \a) have definite energy and momentum (also the spin direction of 
neutrinos is defined) so for the sake of precision we should write 



= \a;p) = a; P,E, a 



(19) 



We can easily construct a state with momentum distributed around a mean value p a . Let 
us assume that, instead of a plane wave, our new state \a) (Eq.(20)) has a Gaussian form, 
which in the momentum representation, is given by 



a-p\a) = m a (p,P a ,cr P p) = 
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(20) 
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where the width a p p is the same for each massive neutrino in the production (P) process 
and the same along all three directions. The average momenta p a of the different mass 
eigenstates are determined by the kinematics of the production process. 

In the wave packet approach, the flavour states \^ a (t)) after time t (given by Eq.(9) in 
the plane wave formalism) are now 



a;P)(a;P\ a 



]T U aa [ d 3 p^ a (p, P a , a pP ) e~ iE ^ 

r, J 



(21) 



where E a {p) = yp 2 + m 2 . 

The same states in the position representation \b; x\ are given by 



X 



b:x) (b: x 



*«»(*)) = J2 U »a[ d 3 x^ a (x, t; v a , a xP 



a; x 



(22) 



where now the function \I/ a (x, t; v* a , a x p) is defined by 

^ a (x, t; v a , a xP ) = J d 3 p^ a (p, P a , a pP ) e~ iEa{p)t 

= ^372 / P> Pa, °PP) ^ 



a;x \ a;P 

X—E a (p)t 



(23) 



Since the Gaussian wave packet in momentum space is picked around the average momentum 
p a we can neglect the spreading of the wave packet and approximate 



E a (p) = E a + V a (p - Pa) 



where 



E a = \fpl + m 2 a and v a ~- 



dE n 



Pa 



Then 



* a (x,t;v a ,a xP ) 
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cxp 



dp P=P« E a 



(p a X -Eat) - V J 



(24) 



>xP 



(25) 



with the width a x p in coordinate space given by 

1 



2a 



P p 



(26) 



As earlier in Eq.(lO), to find the amplitude of the flavour changing process, we project the 
states \ty a (t)} on the flavour states 1^(0)) 



A a ^(t) = (^(0) I tt a (t)) . 



(27) 
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If 1^/3(0)) is the same as before (Eq.(7)) this means that the momentum of each neutrino v a 
is measured precisely. But this is not realistic, so let us assume that also detection process 
is characterized by the spatial coherence width a pD connected with the uncertainties in 
momentum and energy measurements 



l*/?(0)) = Y. U PbJ d3 P^b (P, Pb, <t p p) \b; P) ■ 



(28) 



The average values of the momentum p b are the same as in the incoming wave packets 
Eq.(22). To calculate the spatial decomposition of the detecting flavour state we have to take 
into account that the detector is placed at a distance L from the origin of the coordinates, 
so we have 



= 2 / ^x D \b,x D ){b,x D |^(0)), 



and after the same approximation as before (Eq.(24)) we obtain 



b, x D 



Vp(L)) = Y, U Pb I d 3 x D V b (x* D , 0; v b , a xD ) 

' b J 

(%d, 0; vl, &xd) is given exactly by Eq.(25) after replacements 

a — > b, x^x D , t — > 0, a xP — > a xD . 
The amplitude of the flavour changing process is given by the overlap 



A a ^(L,t) = {*p(L)\* a (t) 



l2a x pcr X D 



UpaUaa exp 



i E a t- p a L 



L — v a t 



(29) 



(30) 



(31) 



(32) 



We have to remember that the origin of coordinates xd and x are not the same, so 

(b, x D \a, x^ = 5 ab 5 {3) (x* D + L - x^j , (33) 

and we have 

A a ^p(L, t) = U Pa U aa J d 3 X^* a (x - L, 0] V a , CT xD ^ ^ a (x , t] V a , (T xP ) 



where the total production and detection width is now 
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(34) 



(35) 



Everything we have done up to now can be applied either to neutrino or neutral boson oscilla- 
tions. The only assumption about the narrow wave packets in the momentum representation 
(Eq.(24)) can be used in both cases. 

Next we calculate the oscillation probability for neutrinos which are relativistic (p 3> m) 
and following Ref. [22] we approximate 



2 2 
777 TTi 

E a = E + ^, Pa ^E-(l-0^, 



(36) 



and 



1 - 



2E 2 ' 



E is the energy determined by kinematics of the production process for a massless neutrino 
and £ is a dimensionless quantity of order unity. We will see (next Chapter) how the rela- 
tivistic approximation for neutrinos causes that the production and detection processes can 
be factorized out and the standard quantum mechanical approach describes the oscillation 
phenomenon properly. 

In all realistic experiments the distance L is a fixed and known quantity, whereas time 
t is not measured. The quantity which we measure is the time integral of the probability. 
Now the time integral can be done, and it is possible to avoid the not properly legitimated 
(in quantum mechanics) replacement x = vt. 



The time integral for 



A a ^(L,t) 



(Eq.(34)) can be done [22] and after normalization 



(j2p P a -,p(x) = lj, instead of Eq.(ll) we have 



X 



Ml, t) 



dt 



\ U Pb\ 2 \ U aaf + 2^2 \U*p b U aa Up b U* ab COS I 2tT— ^ - Lp ab . a/3 
a a>b V ab / 



The oscillation lengths are the same as before (Eq.(12)), namely 

4irE 



T OSC 

ab ~ Ami' 



Am ab = m a 



m 



b- 



(37) 



(38) 



and L c ° h h known as coherence lengths [23] are given by 

4 V / 2^E 2 



- coh 
J ab 



I Am 2 



(39) 
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Comparing Eq.(37) to the usual expression for the neutrino oscillation probability we can 
see that two additional terms appear. 
The second factor 



is equal to unity if a x <C \L°l c \ . This inequality must be satisfied to observe any oscillation. 
The presence of the term (40) which goes to zero for a x > \L^l c \ , reflects the requirement 
which we qualitatively discussed in the previous Chapter: to see the oscillations, the local- 
ization of the source and the detector must be much better determined than the oscillation 
length. 

The first factor 




was predicted long ago [23]. It is connected with the fact that two wave packets each 
with different momentum and energy, have slightly different group velocities. It means that 
after some time the mass eigenstate wave packets no longer overlap and cannot interfere to 
produce oscillations. It is very easy to predict the value of coherence length. If both wave 
packets have width a x p along the direction of propagation and the difference between group 
velocities is \v a — Vb\ = Av then we can expect that after traveling a distance L 

_ 2a xP v a + v b 
L ~ Av 2 (42) 

both wave packets cease to overlap each other. This L is just the coherence length. 

For relativistic neutrinos, using approximations given by Eq.(36), we reproduce L c ° b h from 
Eq.(39), to a factor of a/2- We can expect that the coherence length becomes longer if the 
spreading of the wave packets is taken into account. It is indeed the case as it was proved 
inRef.[21]. 

From Eq.(39) we see that the coherence length L c ° b h is proportional to a x = ^ a xP + a xD 
and not only to o x p as in Eq.(42). It means that precise measurements of momenta of all 
particles appearing in the detection process (which implies small a P D, thus large <j x d) can 
increase the coherence length [22, 24]. This is a wonderful example of quantum mechanics in 
action. A measurement can restore the coherence. Two wave packets having negligible over- 
lap in the detector (thus without detector influence they cannot interfere, and the oscillation 
disappears, a x = & x p), because of precise measurements (<j x d &xp) , may still interfere 
to give rise to oscillations (a x ^> a x p). This feature of quantum mechanics disagrees with 
causality. However, it is not the first time when quantum mechanics is at variance with 
common sense. 

But Eq.(37) also restores some common sense. Measurements of momenta and energies of 
detected particles cannot be too precise if we want to maintain the particle oscillation. As a 
matter of fact, we have a longer and longer coherence length, but on the other hand increasing 
a x makes the position of the detector to be more and more undefined. If o x > \L^l c \ the wave 
packets of neutrinos v a and v b lose coherence (Eq.(40)), the oscillation between them is wash 
away. We see, particularly, that the plane wave approach and oscillations are incompatible. 
For a plane wave a x — > 00 (o p = 0) and the factor (40) teaches us that oscillations disappear. 

There are also approaches where the neutrino oscillation is treated in a manifestly Lorentz 
invariant way [25] . The final answer is exactly the same as was presented up to now, but for 
one difference. In the fully covariant treatment, besides the spatial width a x , also temporal 
width a t should appear. It causes only one change in the oscillation probability formula 
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Eq.(37). Instead of the spatial width a x , a new effective one, a a b appears in the coherence 
length (39) and in the factor (40) 



a 



X 




V a + V b 



(43) 



where v a and vj, are group velocities of wave packets. 

If a t is given by the lifetime of particles which produce neutrino, o t — r [25] (e.g. pions, 
kaons or muons) then the second term in (43) is usually much bigger than the first one 
and the role of factor (40) becomes more important. We should stress however that the 
independent widths of spatial and temporal characteristics of wave packets cause that freely 
propagating particles are not necessary on mass-shell. If we insist to have our particle 
on mass shell (all the time, not for mean values exclusively, which is equivalent with the 
requirement that our particle's state satisfies the equation of motion) only the momentum (or 
only the energy) distribution should be applied. Then energy (or momentum) is distributed 
also but in agreement with the on-shell relation E = ^p 2 + m 2 . Such an approach, which 
we also used in our presentation above is the standard one. 

IV. CURRENT DENSITY APPROACH TO PARTICLE OSCILLATIONS AND 
THE PROBLEMS OF STATES WITH UNDEFINED MASSES. 

A typical experiment which tries to observe particle oscillations measures the flux of 
flavour /3 type particles in the detector localized at some distance L from the source which 
produces particles with flavour a. The time of measurements is not known. Usually typical 
measurements last hours, days or even years (like the observation of solar neutrinos). So the 
most appropriate way to find the probability (or number of particles) to cross the surface 
dA of the detector (see Fig. (1)) is to integrate the probability current density over the 
surface and integrate the result once more over the duration of measurements 



This procedure seems to be so easy and natural, that we can ask why people try to use 
other, more complicated methods. The answer is very simple, there is a problem with the 
correct definition of the current jp(x,t). This current should be defined for particles which 
we measure, that is K°,l( ,i' e or v^. But these particles have undefined masses, and we 
do not know how to define the probability current for such particles. The problem is more 
general: how to define properly the creation and annihilation operators for undefined mass 
states [26,27]? Here we will not discuss all trials [26] to resolve this problem in quantum 
field theory. We will concentrate only on a simple example of definition of the current jp 
for particles with flavour (3. For relativistic currents the problem has not been solved. For 
kaons it was done in Ref.[18]. In the non- relativistic case the free Schrodinger equation for 
particles with mass m a can be written in the form 




(44) 




(45) 
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which is the appropriate non-relativistic limit of either Klein-Gordon or Dirac equations 

2 

(E ~ §^+ m )- For the Schrodinger equation (Eq.(45)) we know how to define the probability 
current: 



j a (x,t) = — Im(V* a (x,t) V V a (x,t)) 

for which the usual continuity equation is satisfied 

d 



dt 



+ div j a = 0. 



The problem arises while we try to define the current for the states 



cm a + dm b and Vp = -d*^ a + c^! b , \c\ 2 + \d\ 2 = 1 



(46) 



(47) 



(48) 



which are an orthogonal mixture of two states with different masses m a and m b . We expect 
that, because of mass mixing, the currents for ^ a ,p will not be conserved [18], so let us 
propose a modified "continuity equation" for ^ a ,p states in the form 



d_ 

dt 



I* 



+ div j a: {3= d at p. 



(49) 



With the following requirements, concerning the new current ^^[18], that: (1) only "ve- 
locity" terms with one gradient are included, (2) for m a — > m b . the "diffusion" terms d a ^ 

should vanish d a: p — * , and (3) the sum of both flavour currents j a + jp is conserved 

d 



dt 



(\V a \ 2 + \Vp\ 2 ) + div(j a + jp) = 0, 



the currents j a> p and the diffusion terms d a> p can be found. 



ja= |c| 2 ja + \d\ 2 jb +Im 



1 1 

cd*{ — %grad^ a ^* a grad%) 

m a m b 



and 



d a = (m a - m b ) Im 



cd*{2^ a m b 



m a m b 



■{grad^ a ){grad%)) 



(50) 



(51) 



(52) 



For jp and dp we have 

jp=3a ( c ~d*, d^ c*) , and dp = d a (c^ -d*, d -> c*) . (53) 

Calculations for — K° mixing probability, using the definition (Eq.(44)) have been done 
in Ref.[18]. The standard formula (Eq.(6)), with the same oscillation frequency uo = ^ was 
recovered, supporting previous results. 

There is, however, one objection concerning this approach. The flavour currents which 
we use are not conserved. This nonconservation is given by the diffusion term (Eq.(52)) 
which is proportional to Am ab = m a — m b . So, we can expect that all our calculations 
have been done with the some precision. Then, if the result is proportional to Am ab , our 
probability flux calculations give the correct answer. The diffusion terms will change the 
result in higher powers of Am ab . 
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V. TREATMENT OF PARTICLE OSCILLATION IN THE FRAMEWORK OF 

QUANTUM FIELD THEORY. 



There are several papers where authors investigate the neutrino oscillation problem in 
the framework of quantum field theory [19,20,24,28,29,30]. What are the main reasons for 
those studies? 

First until now, to describe particle oscillations, we have used states with undefined 
masses (Eq.(4) for kaons and Eq.(7) for neutrinos). But, as we have seen in the previous 
Chapter, there is a problem of proper definition of such states [26]. Only in the extremely 
relativistic limit the flavour states are defined correctly [27]. As we will see, in quantum 
field theory the particle oscillation can be treated without resort to weak eigenstates. 

Secondly, we have completely neglected the effect of the production and detection pro- 
cesses. It has been shown [27], that the neutrino oscillation probability is independent from 
the details of the production and detection processes only in the case of extremely relativistic 
neutrinos. 

And finally, in real neutrino oscillation experiments only associated particles, hadrons 
and charged leptons are observed. Neutrinos are not prepared and not observed directly. 
One can measure the energy and momentum distributions of other particles which appear in 
the production and detection processes. As we will see, only the quantum field theoretical 
approach gives the opportunity to express the neutrino oscillation probability in terms of 
measured quantities. 

Let us now describe briefly how the particle oscillation is treated in field theory. As in 
Ref. [20,30] we will describe the process of neutrino production (P) and neutrino detection 
(D) as one Feynman diagram with a virtual neutrino propagating itself on macroscopic 
distances between the source and the detector. Let us consider the process [20,30] 

Pi - Pf + C + "a 1 - ' up + D T - D F + l~ p (54) 

where Pj and P F (-D/ and D F ) are the particles in the production (detection) processes. 

The production and detection processes are localized in coordinates x P (x D ) and times 
tp(tr>) (see Fig. 2). All initial (P/,Pp,Z+) and final (^Dj, D F ,lp^J particles are described by 
wave packets. Their shapes depend on the measurement precision in the production and 
detection processes. The amplitudes for the full process can be written in the form 

Ap^ d = (P F , Z+; D F , IplSlPjDj) . (55) 

We can see that there are no neutrinos in the initial and final states. Only particles which 
really appear in the production and detection "equipments" are observed. Neutrinos with 
mass m a propagate virtually between the source and the detector and are described by the 
Feynman propagators 

(0|T („. („) Vi W ) |0) = / l^^ii^-e-C.-). (56) 

We will not present the details of all calculations, which are straightforward but tedious. 
A clear presentation can be found in Refs. [20, 28, 30]. We will concentrate only on the 
discussion of the final results. 
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First of all, neutrinos are not directly present in Eq.(55), but this is not necessary. The 
amplitude Ap^p, depends on points (xp,tp^ and (xp,,tp,^j where neutrinos were born and 
detected, and this is enough to study oscillations. 

Next, the amplitude Ap^p, depends also on amplitudes of the production and detection 
processes and the full structure of A P ^ D is the following 

Ap^ D = U Pa U aa A a f a [x* D - Xp, t D - t P ) , (57) 

a 

where A a describes the process of neutrino creation and annihilation. The standard oscil- 
lation formula is recovered only if A a can be factorized. This happens, when amplitudes 
A a become independent of neutrino masses, A a = A . If all neutrinos are relativistic then 
A a = A (m a = 0) and the oscillation probability can be defined. 

In case of relativistic neutrinos the time integrated neutrino flavour changing probability 
is given by the similar formula to Eq.(37) with two changes. First, the dumping factor 
(Eq.(40)) is slightly modified and now equals 

e ^ L ab J (58) 

where £ is the same quantity as before (Eq.(36)), but u is the new factor which depends 
on the production and detection dynamics and can be large (e.g. u = 10 is possible [30]). 
The second modification is a little different definition of the coherence length. Instead of 
(Eq.(39)) we now have 

4E 2 - 

'aft I 



with the same factor u as in Eq.(58). 

The oscillation length L ^ b c and the spatial width a x are given by the same formulae as 
before (Eqs.(38) and (35)). 

There is also an additional very important difference between the present and the former 
wave packet approaches. Before, a x p and o x p> (Eq.(35) were two spatial widths of neutrinos 
specified in some way by the production and detection processes respectively. Now, these 
quantities are defined by spatial widths of hadrons and leptons which are measured. It turns 
out (for detail see Ref.[30]) 

1111 



a xP (T \p 1 a lp F a xa 



and 



+ + (6°) 



a xD (j2 xD l a lp F a x(3 



The widths of observed particles in the production and detection processes define the width 
of " the neutrino" , even if there is no place in this approach for physical neutrinos (only 
virtual ones appear). In the configuration space (Eq.(60)) the sum of the inverse squares 
of widths for the observed particles gives the inverse square of the resultant width. Then 
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the smallest ingredient width dominates the values of a x p or <t x d- It is the opposite for the 
resultant width a x , where a 2 = a 2 P + u 2 D . From the definition of the momentum width 
Up = ^- , it follows that it is just opposite in momentum space, so then 

1 1 1 

2 2 2 ' 

<j p a pP a pD 

but 

°Ip = °l Pl + °Ip f + (62) 

and the same for particles in the detection process. We can see explicitly, that precise 
measurements of momenta of all particles involved in the neutrino detection process (small 
VpDn&pDp and <7 pa ) give a small resultant width <j p d, thus large <j x d and large o x . The 
same subtle thing, which we have discussed before in Chapter 3, that the final measurement 
is able to recover the interference, in the present interpretation has found a much stronger 
background. 



(61) 



VI. DO PARTICLES RECOILING AGAINST MIXED STATES OSCILLATE? 

For many years oscillations of particles like kaons or neutrinos were treated in isolation. 
The circumstances in which oscillating particles were produced have not been considered. 
Recently a series of papers has appeared [31, 32], in which the kinematics of the production 
process have been taken into account in detail. Authors claim that, because the produced 
oscillating particles have neither momentum nor energy defined in explicit way, this fact 
should have consequences not only for them but also for the recoiling particles. 

Let us consider the K° production in the reaction 

n- p -> AK°, (63) 
or the neutrino production in the tt + decay 

7T+ -> fJ+Vp. (64) 

If the invariant mass of the initial system is denoted by M (M 2 = (P w - + P p ) 2 or M 2 = m 2 ) 
then energies {Ej) and momenta (p^) of outgoing particles depend on masses (mi) of K° or 
v^. In the CM system there is 



Pi 



1 /2 

(M 2 — m 2 — m 2 ) — 4m 2 m 2 



2M 
and 

Ei = M2 + m *~ m \ (65) 

2M v ; 

where m is the mass of the recoiling particle A or From the four-momentum conservation 

in the production processes (63) or (64), the energy and momentum of A or are also 
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defined by (65). Authors of Ref. [31] claim that if K° or oscillate, also the recoiling 
particles A or fi + do the same. If it is true, the existence of such phenomena could give a 
chance for indirect observation of neutrino oscillations which are very difficult to observe in 
a direct way. However other papers have immediately appeared [33, 34, 35] where authors 
have been strongly against the oscillation of particles produced in association with kaons or 
neutrinos. We will present here our approach to the problem [36] which also supports the 
opinion against a visible oscillation of the associated particles. To fix notations, everything 
will be described for the process (63), but equally well we can show the lack of visible muon 
oscillations in the pion decay (64). 

First of all we would like to specify the kind of oscillation, we can consider for A (or 
fi + ). Two A's with different masses do not exist. But even without mass differences the 
A's are produced in association with the long-live Kl and the short-live Ks- As Kl and Kg 
have different masses A's will be produced in two orthogonal states with different energy 
and momentum. 



A L ) = -p L ,M- E L 



and 



|Ac) = 



p s , M-E s 



(66) 



where J? us) an< ^ El(s) are the momentum and the energy of the K L (K S ) in total CM frame 
of the process (63). We do not know in which state \A L ) or \A S ) the A particles are produced, 
so let us assume that at t — they are produced in some state which is a linear combination 
of both states (66) 



|A(0)> =a|A L ) + 6|A s ) 



I I 2 , 

a + 



(67) 



As both ingredient states have different energy they evolve with time in a different way, so 
there is some chance that after some period of time the state |A (0)) will oscillate to the 
orthogonal one 



A' (0) 



-b* \A L ) + a* \As) 



(68) 



Do we have a chance to recognize both states |A (0)) and A' (0)y? In the neutral bosons 

system, because of the strangeness conservation, K° and A°interact strongly in a completely 
different way and are easily distinguishable. Here we have the same particle A with only one 
decay width I\. In spite of that, in principle we can distinguish |A (0)) from A' (0)y but 
in a much more sophisticated way. A's in both states will decay in the same way (mostly 
to pn~). But because two states Al(s)) have slightly different momenta (in CM frame) also 

will be slightly different. As two states |A (0)) 
, the angular distribution of the protons in CM 



the angular distribution (e.g. for protons) 



and 



A' (0)) are various mixtures of 



A 



L(S) 



frame which come from |A (0)) or from A' (0)j will be different. 

We can see that in principle both states (67) and (68) could be distinguished. But do we 
have anything which may be distinguished? In other words, if we have flux of A's produced 
with kaons in the reaction (63), will their number in the state |A (0)) or A (0)) change with 



distance from the reaction point? For simplicity we will present the answer to this question 
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in the plane wave language. We know that it is not precise, but the value of oscillation 
length obtained in this way is correct. The full wave packet approach, together with particle 
correlations (the EPR effect) will be presented in the next Section. 

Let us assume that at t — 0, A's are produced in the pure state |A (0)) (in the reaction 
(63), the coefficients a = b = ^= , but it is more transparent to leave them undefined). 
We will consider the production of A's and the decay A — > pn~ together. The amplitude 
for A production and decay after a period of time t can be written in a form, where two 
indistinguishable ways of reaching the final state are added coherently [37] 



A (A -> prr") = A (A -> A L ) e -< mA -^) TL A (A L -> pn 

+ A (A -> A s ) e^K-i^M (a s -> jwr") , (69) 



where A (A — > A Lj s) are amplitudes for A production in the states A L ( S )), A (A L:S — > £>7r~) 
are decay amplitudes from both states, m r and T are mass and decay width of the A. The 
different proper times which elapse in the A's rest frames during the propagation are the 
crucial points in our discussion. 
If we denote 

A }fl S ^ W \ = Vsl = \vsl\ e ip ^,a = \a\ e*«, b = \b\ j», (70) 
A (A L -> pn ) 

the probability for A production and decay from the initial state |A (0)) can be written 

P (A -> = \A (A -> p7r~) I 2 = U (A L -> pvr") I 2 (71) 

x {| a | 2 e- r ^ + |&| 2 |775L| 2 e- rT - s 

+ 2 |a% 5 | 2( T i+ r s)r cog ^ _ T ^ + (fb + pLS _ (/?a ]} 

The oscillation can possibly arise from the term m\ (t l — t s ). If instead of A we consider the 
production and decay of the initial kaon a similar formula would be obtained but with one, 
as we will see, crucial difference. As masses of K L — K s bosons are different, the oscillation 
factor is equal to ttilTl — msTs- 

How to calculate the proper times? They are measured in different Lorentz frames, in 
the rest frames of A L and A s . The basic principle of quantum mechanics - the superposition 
principle tells, that we can add two states at the same time 

|*(0> = M0> + M*)> ■ (72) 

In the position representation we add wave functions 

(x\*(t)) = {x\<p 1 (t)) + {x\<p 2 (t)) (73) 

in the same position x and at the same time t. It means that the proper times are not 
suitable variables. We have to add wave function at the same point (x, tj of the same 
Lorentz system. It is necessary to transform tl and t$ to the same common frame. The 
CM frame for the whole 7i~p — > A°K° process is the most convenient in this place. 
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For convenience, we consider only "the one dimensional" problem x— (a;, 0,0). Then 
Lorentz transformations between the rest frames for Al(s) and the CM frame are given by 

tl = lh(t- (5 L x) , £ L = 7 L (x - p L t) , 

and 

rs = 75 (t - Psx) , £s = Is (x - Pst) , 

where 

W = f 2 . ft* = fL. (74) 

At the beginning t — 0,x — 0, and two "ingredients" of the A particle, and As are created 
at tl = rs = and £l = £s = 0. But particles in two states A^ and As have different speeds 
and after time t they are in different points in the CM frame (see Fig. 3). 

In classical mechanics, for point particles, it is impossible to have a situation that two 
particles which were born in the same point and at the same time but moving with different 
speeds would be still in the same, common points at the same time later (Fig.3). Accordingly 
to our previous statement (Eq.(73)) such particles will not interfere for any time t > 0. But 
in QM particles are described by wave packets (in the limiting case-plane waves). We do not 
know at what place the particle was born inside the wave packet and what was the speed of 
it (see Fig. 4). It is not strange that different parts of two wave packets still interfere. Inside 
wave packets, energy and momentum are distributed in agreement with QM prescriptions 
and it is not a surprise that they can be not conserved [32]. 

As we remember, in the wave packet approach, to find the probability as a function of 
position it is not necessary to assume any relation between t and x. We simply integrate 
over t [36]. It is possible, however, to find such a relation between t and x that the oscillation 
length, which we obtain in this frame, will be (to the first order) the same as in a proper 
wave packet approach. Such a frame was found [35]. It is the CM frame for A^ and As, 
where their momenta are opposite p* L = — p* s (see Fig. 4). The velocity (3 of the origin of 
Al — As center of mass frame in the laboratory system can be easily found from the relation 

7 (p L - PE L ) =p* l = -p s = -7 (PS - PEs) , (75) 

so 

P= PL+PS = El ~ E s. (76) 
El + Es Pl-Ps 

Then the movement of the origin of the A^ — A^ CM frame in our laboratory system is 
described by the obvious relation 

x = pt (77) 

with P given by Eq.(76). 

Using this classical relation Eq.(77), the proper times tl and t s (from Eq.(74)) are given 

by 
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r L = 1l {j ~ Pl) x, 

(78) 

Ts = 75 (| - Ps) x 



and taking Eqs. (74) and (77) we have [35] 



t l -t s = x 



1l\^-Pl\-1s (j}-Pl 



(79) 



and the oscillation length (Eq.(71)) is infinitely large. If we calculate, in a analogous frame, 
the oscillation factor for kaons we obtain [35] 

m L T L - m S T S = m L -f KL I—-/3 Kl \ - m s ^ Ks ll--/3 Ks j 



m 2 L -m 2 s _ mL + ms (m L - m s ) _ mAm _ Am 



2 



Pk l +Pk s (pk l +Pk s )/2 P 2p 

(80) 

which reproduces the well known result for the oscillation frequency (Eq.(6)) 

x mAm Am 2 

Amt = Am— = x = x. (81) 

v mv 2p 

The relations (79) and (80) which characterize the oscillation length are obtained in this 
special Lorentz frame. In other frames, these results are correct only to first order in Am. 
Instead of checking the dependence of the oscillation length on the Lorentz frame, in the next 
Section we will present the more complete wave packet approach. Here we have found that 
the oscillation length of particles recoiling against mixed states is very large. It means that 
even if we consider the oscillation of such particles (A or //) separately, without connection to 
kaons or neutrinos it is impossible to observe such oscillations on any acceptable terrestrial 
distance. Now we consider the oscillation of both particles (A and K or \x and v) together. 



VII. CORRELATIONS FOR TWO OSCILLATING PARTICLES, EPR EFFECT. 

Up to now we have considered the oscillation of one particle without taking into account 
possible correlations which may appear for two or more particles in the final states from which 
at least one oscillates in the traditional way. There are many such cases. Some of them, with 
one oscillating particle, have been discussed above (n~p — > A°fT° or n + — > /i + f At ) . There are 
also interesting processes with two oscillating bosons e.g. $ — > K°K , \1> (3770) — > D°D or 
e + e- -> T (4s) -> B°B°. 

Let us describe the last of them. At t — the state of two bosons is a combination of 
states with definite masses 

B °B°) t=0 = VabRBaR-Bb \ B aB b ) t=Q , (82) 
a,b 
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where Rbo, and R-^ b are elements of unitary matrix which describes the mixing. The momen- 
tum conservation in the production process gives altogether w(w 2 +1 - > independent momenta for 
all n 2 pairs B a B b . Sometimes there are additional correlations between various mass states 
B a and B b in Eq.(82). If, for example, the B°B pairs are produced by the T(4s) decay 
then the state 



B°B°^ must be totally antisymmetric [38,39] (since T(4s) has intrinsic spin 
s = 1 but B mesons are spinless, the B pair is in a p wave). The factors r] ab in Eq.(71) are 
responsible for such correlations (see for details Ref.[36]). 

Each state 
is given by 



B a {p)) is described by a wave packet which in the momentum representation 



a,P), 



[a,P \B a ) t=Q = J d 3 p^ a (p,P a ,a. 
where fy a (p, P a , a p ) for simplicity is taken as the Gauss function 



(83) 



(P, Pa, (Tp 



2na r , 



3/2 



(84) 



with p a - the average momentum and a p - the width of the distribution. After time t (taking 
into account the particle decay) the state \B a ) will evolve into 



\B, 



a/t=0 



\B a (t))^e- iHt 



B a) t=C 



where r a = J- 



Ea \ Pa 



= e~ik J d 3 p^ a (p, P a , a p ) e- iE ^ \a,p), 



is the lifetime of the "a" particle in a chosen Lorentz frame. 



(85) 



The states \B a (t)) in the position representation will be given by 



\B a it)) = e 2*a J d 3 x^ a (x, t; t£, a BP 



a,x), 



(86) 



where \l/ a (^x,t;v^, a bp) is the Fourier transform of the momentum distribution (Eq.(84)) 

and is given by Eq.(25). Let us assume that two detectors are placed at points L\ and 

L 2 . The detectors will measure the particles with beauty "1" and "2" (1,2 = B°,B°, 
respectively). 

The states of the B mesons measured by the two detectors are defined by 



B 1 (L^j ^ = Y, R icJ d 3 Xl m c (V x - L u 0; v c , a 1D ^j 



C,Xi 



and 



Bo L 



(87) 



(88) 
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Notations in Eqs. (87) and (88) are similar to previously presented in Eq.(30). 

We can find the amplitude for two oscillating particles in the same way as before (Chapter 
3). Then the amplitude of the process where two particles B° and B° produced at t — at 
point x— are detected as particles with beauty "1" ("2") at point Li(L 2 ) at t — ts (%), 
is the following 

^ Bl ^ B2 (4* B ;4%) = (bi (li),b 2 (l 2 ) I B°(t B ),B°(t^ 

= Nj2R*laRBaR* 2b RBbV ab e~^ + ^ 
a,b 



(89) 



with the normalization factor 



N = 



^1D^2DCTBP^ bP 



nl/2 



(90) 



and the effective total widths 



<?B 



a 



ID 



a 



BPi 



a. 



2D 



a- 



BP 



(91) 



The amplitude (89) can be used in various situations. If we apply the formula (89) to the 
description of the EPR effect in the T (4s) — > BB decay (Refs [38,39]), the collapse of 
the BB wave function is included in a natural way. Our approach is an alternative to the 
amplitude description (Refs [38,39]) and, in case of particle mixing, takes into account the 
EPR correlations in a much more transparent way (see Ref.[36] for detail). 

The formula (89) can also be used for the /xz/ M pair "oscillation" from the 7r — > fiu^ decay. 
In this case only one particle, the neutrino, mixes. Then, in the application of Eq.(89) to 
our present purpose, we have to take one diagonal mixing matrix (e.g. R 2 b = S 2 b, R^b = ^m) 
and r/ab = 5 a b- Usually neutrinos are considered as stable or very long living particles. Let 
us also assume that the "oscillation" of the muons is measured on a distance much shorter 
than their decay length. Then both factors in Eq.(89), which are responsible for particle 
decay, may be neglected. In such circumstances, the probability that neutrinos produced as 

= a type together with muons at x= are observed as a /3-type neutrino at distance 
L u and the muons at distance L^, 
Ref.[36]), 



after integrating over times is given by (for details see 



Ppa (L^, L u ) — 



E 



E 

ab 



-2ni 



j^(j,coh 
ab 



V 2 ■ 
u [ia 



AT", 

ab 



-2ni 



x e 



( L » ) 



ab 



(92) 
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The oscillation (L^ sc )and coherence [L^l° h ^j lengths, and the factor N^ b for muons are (for 
neutrinos the appropriate expressions are similar) 



(E^ - E^ b ) Vt T ^ V f - (v^a - Pfib) 

V fia "+" V ph 



(93) 



/ 7 ,2 _i_ 7 ,2 \V2 

LT = 2(7^ ^1^L_ , (94) 

V Ka - V,) / 



and 



- ^) 2 + T&) 



l2' 



(95) 



(£^ a - £^) + tj m6 ) - (p Ma - P/lb ) (tj^ + vfy 

where E^ a ^,p^ b ) and t^ a (b) are energy, momentum and velocity of the muon associated 
with the neutrino a(b). 

First of all, we can see from Eq.(92) that the muon oscillation disappears if we do not 
measure separately the neutrinos with flavour (3. In such case the probability given by 
amplitude (92) is constant in L M and L u and can be normalized 

1 £P Pa (L li ,L l/ ) = l. (96) 



But even if we measure the /3-type neutrino, the muon oscillation will not be seen. We can 
prove this statement because we know precisely the muon oscillation length (Eq.(93)). 
If we denote the difference between the masses of two neutrinos a and b, as 

m a -m b = Am ab , (97) 

the inverse of the oscillation length may be decomposed in powers of Am ab . For muons, the 
term proportional to the first power of Am ab vanishes, and 

(2tt) (LITT 1 = -2C 2 (l - vl) Va 1 (Am ab ) 2 + (98) 
where ( is the factor in the decomposition 

W = *(l+c(^)+/>(^V + ...|, (99) 
while for neutrinos there is 



lab 
Pa 



A™2 

(2vr) (L—r 1 = + •- (ioo) 

and we reconstruct the previous formula for mixing particle oscillation length Eqs.(12), (38). 

From Eq.(98) it follows that L^ sc is very large and for the acceptable neutrino mass 
difference the muon oscillation length is much bigger than its decay length 

L T C > ct m ~ 660m. (101) 

We can see that even if the neutrino and the muon are both measured, the oscillation of 
muon will not be observed. Taking into account opinions presented in the latest exchange of 
views [31-35] we agree with the statement that, in practice, oscillations of /i or A particle 
are impossible to observe. 
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VIII. SUMMARY 



Let us briefly summarize the main results of this review paper 

(i) Many conceptual difficulties arise in the plane wave approach to the particle oscillation 
problem. This approach gives us the shortest way to get the correct expression for the 
oscillation length, but it fails if we try to describe other aspects of the particle oscillation. 

(ii) In real oscillation experiments neither energy nor momentum are the same for all 
eigenmass particles. 

(iii) The wave packet approach 

- gives the proper oscillation length L°l°, 

- introduces the concept of the coherence length L^f 1 , such that for distances greater 
than LrfJ 1 the particle oscillation disappears, 

- in the proper way, takes into account the fact that to observe oscillation, the sizes of 
particles source and detector must be much smaller than the oscillation length, 

- gives a possibility to understand in a simple way the phenomenon, that a precise 
measurement of the detected particles momenta may restore the coherence between various 
eigenmass states and, as a consequence, the oscillation between particles, 

- temporal and spatial distributions in the wave packet are correlated by the requirement, 
that particles are on mass shell. Independent distributions for time and space give a wave 
packet which does not satisfy the free particle wave equation. 

(iv) A problem arises with the proper definition of the Fock space for flavour states. 
As a consequence, the defined probability currents for such states are not conserved. Then 
the calculated flavour changing probability is correct only to the first power of the mass 
difference Am. 

(v) The most adequate approach to particle oscillation is given by quantum field theory. 
It can be seen especially for neutrinos, which are "neither prepared nor observed", and only 
propagate between sources and detectors. In this approach 

- production and detection processes are fully taken into account as in real experiments, 

- physical quantities are expressed in terms of measured quantities, like momenta and 
energies of hadrons or charged leptons in the neutrino creation and detection processes, 

- flavour states, which are not well defined, are not necessary to describe the oscillation 
process. Only the mass eigenstates and the elements of the mixing matrices may be used, 

- it is clear in which circumstances the production and detection amplitudes can be 
factorized out and the separate oscillation probability be defined, 

- the oscillation of non- relativistic particles can be described in the proper way. 

(vi) We have specified the meaning of oscillation of particles which recoil against mixed 
states (as A in the process ir~p — > A A" or \x in the decay 7r + — > p^v^). Even if we 
consider oscillations of such particles separately, without connection to kaons or neutrinos, 
it is impossible to observe these oscillations on acceptable terrestrial distances. 

(vii) We have found the general amplitudes which describe the oscillation of two particles 
in the final states. These amplitudes can be applied to 

- description of the EPR correlations in decays like T(4s) — > BB or $ — > KK, including 
the mysterious collapse of the wave function in a natural way and giving the possibility to 
discuss the relativistic EPR correlations on distances longer than coherence lengths. 
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- description of two particles oscillation from which only one has indeterminate mass like 
AK° or fiu^ . Oscillations of particles with known mass (eq. A or //) can be defined only if, 
in the same time, flavours of the unknown mass particles are measured (K° or v^). In this 
case, however, the oscillation length of particles with determinate mass is very large, much 
larger than the particle decay length, which makes it impossible to observe their oscillation 
in practice. 

I thank J. Sladkowski for valuable remarks and M. Czakon and J. Gluza for reading the 
text. 
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FIGURES 




Detector 



FIG. 1. Particles with the flavour a are produced in the "Source". After travelling the distance 
L they are detected as particles with flavour (3. The probability of such detection is given by the 
probability current integral over the detector active surface (6 A) and over the time of measurements. 




FIG. 2. In the field theory approach to the neutrino oscillation only the source i 
and detector (Dj, Dp, e~) particles are measured. Between the production (xp,tp) and detection 
{xD,to) points neutrino propagates as virtual particle. 
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FIG. 3. Relation between CM frame for the production process ir~p — > K°A° and the two 
rest frames for the and As which move with different speeds. In classical mechanics for t > 
the Al and As are in different points. 
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A L + A 5 CM frame 



FIG. 4. The same as in Fig. 3 but in quantum mechanics where the Ax, and Ag particles are 
described by the wave packets. Even if the centres of the wave packets at the same time are in 
the different points (like for classical particle in Fig. 3) the two states for Ax and Ag may interfere. 
The interference will be possible if the two wave packets overlap. 



28 



